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compared to the CP-conserving case. The branching ratios exhibit, however, a strong depen- 
dence on the choice of the phases in the Lagrangian which affects the significance of flavor 
models. In particular, for general CP-phases, the lepton flavor violating rates appear to be 
essentially uncorrelated with the possible high- and low-energy lepton mixing parameters, 
such as the reactor angle. 
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1 Introduction 



Neutrino oscillation experiments have during the past decade pinned down the neutrino mass 
and mixing parameters to a remarkable precision [U El EJ H] . Global fits [5] tell us that the 
solar and atmospheric neutrino mass squared differences are (at la) 

Ami = 7-651^ x l(T 5 eV 2 , |Am 2 t J = 2.40+°;^ x l(T 3 eV 2 , (1) 

whereas the solar, atmospheric and reactoiQ mixing angle are respectively given by 

0i 2 = (32.6± 3 ; 2 )°, 23 = (45.0+^3)°, 13 < 13.2° (3a). (2) 

The leptonic Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing matrix £/pmns [2] can be 
well approximated by the Harrison-Perkins-Scott (HPS) tribimaximal mixing matrix £7hps 
[7J (up to phases) as 



f^PMNS ~ f^HPS 



1 1 1 

V6 V3 V2 

\ J- L _!_/ 



(3) 



In ?7hps, the solar and the atmospheric angle are given by 0i2 ~ 35° and 023 = 45°, whereas 
the reactor angle 0i3 vanishes. The measured PMNS mixing angles may thus be treated 
as deviations from exact tribimaximal mixing [SI E] , describing nearly tribimaximal lepton 
mixing [10J. 

The smallness of the absolute neutrino mass scale m v ~ 5 x 1CT 2 eV can be understood 
in terms of seesaw mechanism [HI [12] , which establishes a connection between low-energy 
neutrino observables and physics near the unification scale Mx ~ 2.5 x 10 16 GeV [13J. In 
such a context, neutrino oscillations can therefore probe high-scale theories of flavor. 

It is interesting to ask whether there are possibilities to test mass and lepton flavor models 
in other than only neutrino oscillation experiments. Such a possibility is offered by charged 
lepton flavor violation (LFV) in supersymmetry (SUSY). In the standard model (SM), LFV 
is absent and adding right-handed neutrinos to the SM only tiny LFV effects that are sup- 
pressed by the smallness of the neutrino masses. In SUSY, however, virtual effects of the 
right-handed neutrinos and their superpartners affect the renormalization group equations 
(RGE) of the slepton masses and trilinear couplings and induce observable branching ratios 
(BRs) of LFV decays [H]. 

One attractive feature of the seesaw mechanism is that it allows for leptogenesis via 
the decay of the right-handed neutrinos [15J (for reviews see [16]). Barring special cases, 
however, the complex Yukawa couplings providing the necessary CP-violating phases can 
generally lead to a large enhancement of the LFV rates by large factors > 10 [17]. If these 
CP phases are not fixed, their presence may obscure possible predictions from flavor models 
and, furthermore, make it difficult to distinguish experimentally different candidate theories 
of flavor. For a discussion of LFV in certain quark-lepton complementarity scenarios see for 
example [IS] . 

1 Note, the current best-fit value for the reactor angle is #13 = 6.5°. 



2 



In this paper, we study the LFV decays /i — » cy, r — >■ 7x7, and r — >■ e7, in minimal 
supergravity (mSUGRA) for a large class of CP-conserving and CP-violating lepton mass 
matrix textures that yield nearly tribimaximal lepton mixing. The textures are formulated 
at the level of the Lagrangian and are characterized by hierarchical entries that can be 
generated in flavor models using, e.g. the Froggatt-Nielsen (FN) mechanism [19]. The tex- 
tures implement the idea of quark-lepton complementarity (QLC) [20] in a general way: 
They realize nearly tribimaximal lepton mixing by taking large contributions from differ- 
ent lepton sectors, i.e. the charged lepton and the left- and right-handed neutrino sector, 
into account. The FN mechanism has, for instance, been used in Refs. [2D 1221 1231 125] to 
construct charged lepton and even neutrino mass textures which can be implemented by 
discrete flavor symmetries [231 EH EE]. Many models have been proposed in the literature 
to reproduce tribimaximal leptonic mixing using non-Abelian flavor symmetries (for early 
models based on At ~ Z 3 x (Z 2 x Z 2 ) and examples using the double covering group of A 4 , 
see Refs. [2Z] and [2H])- We demonstrate that the inclusion of random CP- violating phases 
at the level of the Lagrangian can enhance the branching ratios typically up to 2-3 orders 
of magnitude. Moreover, we find a strong dependence of the LFV-rates on the choice of 
the phases in the Lagrangian. As a consequence, it appears that the assignment of random 
CP-violating phases in the Lagrangian in general erase possible correlations between LFV 
rates and the PMNS mixing parameters of the low scale theory. 

The paper is organized as follows: In Sec. |2j we introduce the SUSY seesaw mechanism 
and its different parametrizations. Charged lepton flavor violation in SUSY is described in 
Sec. [3} The types of textures for which we study LFV are introduced in Sec. |4j and our 
results for LFV are then presented in Sec. |5j Finally, in Sec. [6j we give our summary and 
conclusions. 

2 SUSY Seesaw Mechanism 

In what follows, we assume that the left-handed active neutrinos acquire their masses via the 
type-I seesaw mechanism [UJ. We consider the seesaw mechanism in SUSY, which differs 
from the usual seesaw mechanism in the respect that it involves two Higgs doublets. In the 
type-I seesaw mechanism, the part of the superpotential generating the charged lepton and 
neutrino masses is given by 

W = -(Y^e^H, - (Y D ) t ^H 2 + ^(M R )^, (4) 

where = [y^ ej) T , e°, and z/ c (z = 1, 2, 3 is the generation index) are the matter superfields of 
the left-handed lepton doublets, right-handed charged leptons, and right-handed SM singlet 
neutrinos, respectively. In (Jil), H\ and H 2 are the usual Higgs superfield doublets generating 
the down- and up-type masses, respectively, and Yg and Yd are the 3x3 Dirac Yukawa 
coupling matrices of the charged leptons and neutrinos, whereas Mr is the 3x3 Major ana 
mass matrix of the right-handed neutrinos. After electroweak symmetry breaking, the Higgs 
doublets develop vacuum expectation values (H®), where {H®) = vsm(3 with v = 174 GeV 
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and tan/3 = (i?" ) • The resulting lepton mass terms become 



-(M^-efo - (MD^vfa + ^(Mr)^^ + h.c. 



(5) 



where = (H®)Ye is the charged lepton and = (H%)Yd the Dirac neutrino mass 
matrix. Mi and are complex 3x3 matrices. The matrix Mr is complex and symmetric 
and has matrix elements of the order of the B — L breaking scale Mb-l ~ 10 14 GeV. After 
integrating out the heavy right-handed neutrinos, we obtain the effective low-energy 3x3 
neutrino Majorana mass matrix 



M eff = M T D M^M D = Y£Mr X Y d {v sin (5)\ 



(6) 



giving masses ~ 10~ 2 eV to the light neutrinos in agreement with observation. The seesaw 
mechanism is attractive since it establishes a connection between the absolute neutrino mass 
scale and M B _ L , which is close to the grand unified theory (GUT) scale M x ~ 10 16 GeV. 

The mass terms in (5|6) are diagonalized by unitary 3x3 matrices U x which in their most 
general form can be written as (cf. Refs. j2TJ 122] ) 



U, 



x 



t',D,D',R,p, 



(7a) 



where D x and K x are given by D x = diag(e 1 ^ 1 , e 1(p2 , e 1 ^ 3 ) and K x = diag(e 1Ql , e ia ' 2 , 1) with 
(fi G [0, 2tt[ and a« G [0,7r[ and 



Cl2Cl3 



— S12C23 — C12S23S13& 
AS 



id 



S12C13 $13 e 
C12C23 — Sl2S23 s 13 el<5 S23C13 



(7b) 



\ S12S23 — Ci2C23Si3e — C12S23 — Sl2C23 s 13 e C23C13 J 



sin 6{j , C{j 



is a unitary matrix in the standard parametrization with 
% G {^12, #13, #23} G [0, |1 , and 5 G [0, 2tt[. The mass matrices then read 



cos 6ij, where 



M t = UpMf * g Uf , M D = U* D ,M^U T D , M R = U R M A R ^Ul, M eS = U u M^Uj, (8) 
where the diagonal mass matrices 



Mf ag 



= diag(m e , m M , m 
diag(mf , m 2 , m 



f), 



diag 



= diag(m x , m 2 , m 3 
M c g ag = diag(mi, m 2 , m 3 ), 



having positive mass eigenvalues. The PMNS mixing matrix is given by 



PMNS 



U e U v — ?7pMNS-^Maj, 



(9a) 
(9b) 

(10) 



where the matrix ?7pmns i n (10) is described in the standard parametrization (7b) by the 
solar angle Qyi = On, the reactor angle #13 = 6*13, the atmospheric angle #23 = #23 and the 
Dirac CP-phase 5 = 5. i^Maj = diag(e 1?il , e 1 ^ 2 , 1) contains the Majorana phases <pi and 02 
with 1>2 G [0,n[ m 



2 The phase matrix has been absorbed into the right-handed charged lepton sector. 
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Rotating the lepton doublets t L by Ue to the basis where Mi is diagonal, M e g is diago- 
nalized by the PMNS matrix 



^PMNS^eff^PMNS ~ M eS & - (H) 

Unless stated otherwise, we will from now on work in the basis where the charged lepton and 
heavy Majorana neutrino mass matrix are diagonal, i.e. where M e = Aff ag and M R = M R ag . 

The significance of the mixing matrices U x in ^ is that the mass terms in ^ may 
be predicted by some flavor model such as the Froggatt-Nielsen mechanism [19] or discrete 
flavor symmetries [231 GH 126] . Therefore, by tracing {/pmns back to the matrices U x , which 
describe the high-energy lepton mixing, one can gain better understanding of the observable 
PMNS parameters in terms of a possible fundamental theory of flavor. 

Flavor models do not fix the absolute scale of the Yukawa coupling matrix. In order to 
properly normalize all couplings, we shall express Yd as a function of the absolute neutrino 
mass scale m v ~ 10~ 2 eV and Mb-l ~ 10 14 GeV. For this purpose, consider arbitrary M e g 
and M R and let 7713 and mf be the heaviest mass eigenvalues of M e g and Mr, respectively. 
We normalize M e g by redefining 

Tfl 

M eS -> M' eS = — M eff . (12) 
m 3 

Adjusted in this way, the heaviest mass eigenvalue of M e g becomes equal to m v . Similarly, we 
re-scale = Mb-l to normalize Mr. The normalization of M e g in (12) is then equivalent 
to the redefinitions 



M r ^M' r = Mb ^M r and M D ^M' D = J m ^ MB j±M D = (Hl)Y D) (13) 



where Yd is now properly normalized. From ([8]) and (7a) we find that in this basis 

M D = K* R U R DU* D , M^KUlDnU; . ( 14) 

For simplicity, we will make this rescaling in the following without noting the prime in M R 



and Mjj. From (14) and (13), we see that 



YD = W)\ ^^rJ/^frryM^mDuU;, (15) 

where we have introduced K = KdK* d , and D = D* D ,D* R . The neutrino Yukawa coupling 
matrix in (15) can also be written as follows [29 Pj 



Yd = — - V< aS " R ■ \Mt S " ^pmns- (16) 

(J olll L) 



5 Note that the definition of ZJpmns in differs from our definition by complex conjugation (cf. (11)) 
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Here, R denotes a complex orthogonal matrix which may be parametrized in terms of 3 
complex angles 9i = Xi + iyt as 



R 




— C1S3 — S1S2C3 S1S3 — C1S2C3 
C1C3 — S1S2S3 — S1C3 — C1S2S3 . , 
sic 2 cic 2 



(17) 



with (q, Si) = (cos 9i, sin 9i) = (cos Xj cosh y$ — i sin Xj sinh?/;, sin x, cosh y^ + i cos Xj sinh?/j). 
The parameters can take the values Xj G [0, 27r[ and ?/j G ] — 00, oof (in practical cases, 
however, the yi are constrained by perturbativity to values \yi\ < 0(1), see Sec. 4.2). While 
the light neutrino masses rrii and the mixing angles 9ij have been measured or constrained, 
the phases 0; and S, the heavy neutrino masses mf and the matrix R are presently unknown. 



From (15) and (16) we obtain 



R = (Mf s )-V^*^^,Mg as ^^^C/;^ aj (M5 ag ) 



-1/2 



(18) 



where we have used (10). The parametrization in (16) has the advantage that (i) one can 
understand the impact of nonzero CP-violating phases on the LFV rates in a compara- 
tively simple way and (ii) one can quickly scan the parameter space while ensuring a valid 
low-energy phenomenology (lepton masses and PMNS angles). We will make use of these 
properties in subsequent sections. Note, however, that the drawback of the parametrization 



(16) is that the exact connection with the lepton mass terms in the Lagrangian has been 



lost after rotating to the basis where Mi and Mr are diagonal. 



3 Charged Lepton Flavor Violation in SUSY 

The heavy neutrino mass eigenstates v\ introduced in the seesaw mechanism are too heavy to 
be observed directly but they influence the mixing of the sleptons in the MSSM via radiative 
corrections. The 6x6 slepton mass matrix may be written as a sum of two parts, 



m ? / ml m%\ + f 5m 2 L &™lr \ (10) 
' V m 2 LR m 2 R J MSSM V 5 m l R 1 



where the first part denoted by MSSM is the usual mass matrix in the MSSM without right- 
handed neutrinos. In our analysis of SUSY LFV processes we adopt the mSUGRA scheme 
of SUSY-breaking, in which case the slepton mass matrix does not contain flavor mixing 
terms. In leading logarithmic approximation the corrections to the slepton mass matrix due 



to right-handed neutrinos, denoted by v c in (19), can then be written as [M 



8ml = -^{Zml + Al)YlLY D , 
3 

16^ 2 



Sm 2 LR = -—A v cos WYlLYn, (20) 



where Ljj = \n{Mx/mf)8ij , mf being the heavy neutrino masses, and m and Aq are the 
universal scalar mass and trilinear coupling, respectively, at Mx- With the neutrino Yukawa 
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Br(/x — > ej) 


Br(r /ry) 


Br(r —7- ej) 


Present 


1.2 x 10- n [30] 


6.8 x 10" 8 [31] 


1.1 x 10" 7 [32] 


Expected 


10" 13 [33] 


« 10" 8 


« 10" 8 



Table 1: Current bounds and expected future sensitivities of direct experimental LFV 
searches. 



matrix Yjj and the heavy neutrino masses mf in a given scenario as input at the GUT scale 



we calculate the slepton mass matrix at the electroweak scale using (20). The flavor off- 



diagonal virtual effects in (20) induced by the mixing in the neutrino sector lead to charged 



LFV. More details on this mechanism can be found in 1 35 1 and the references therein. 



3.1 LFV Rare Decays 

In the SUSY seesaw model considered here, LFV processes mainly occur via intermediate 
left-handed slepton flavor transitions. The most important low-energy processes are the 
rare decays U — > Ijj, i 7^ j G e, fi, r, which provide the most stringent bounds on LFV in 
the SUSY seesaw model as of now. The current bounds on these processes as well as the 
expected sensitivities of future experiments are listed in Table [Tj 

Each LFV transition is suppressed in a given process by a small factor \ {5m 2 L )ij /rh?\ 2 (i 7^ 
j), where (Sm 2 L )ij are the off-diagonal elements of the left-handed slepton mass matrix m\ 



specified in (20 ) and m is of the order of the relevant sparticle masses in the loops involved in 



the process. To lowest order in the off-diagonal mass corrections one has approximately 

Br{k -+ L 7 ) oc a 3 mf \^M. tan 2 ^ (2 1) 

This expression is just used for illustration. In our numerical calculations we use the full one 
loop result for Br(Zj — > Ijj), as given in [35] . 

3.2 LFV Processes at the LHC 

At the LHC, a feasible test of LFV is provided by the production of squarks and gluinos, 
followed by cascade decays of squarks and gluinos via neutralinos and sleptons [36, 37J: 



pp 


— >■ 


Qaqb,gq a ,gg, 


Qa(g) 


-> 


x°2<ia(g), 




-> 




lc 


— >■ 





(22) 

where a, b, c run over all relevant sparticle mass eigenstates. Lepton flavor violation can 
occur in the decay of the second lightest neutralino or the slepton, resulting in different 
lepton flavors, i 7^ j. The total cross section for the signature if If + X can then be written 
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as 



a(pp^X° 2 + X^l?lJxl + X) 
2a(pp -> gg)Bx{g -> 

a 

+ a (PP [ Br (& 9X2) + Br(% -> gxS)] f x ^(xl ^Jx?), (23) 

0,6 J 

with Br(g — >■ = 12 a ^ r (9 ~~ * QQa)Br(q a — > qxQ) an d can involve jets, leptons and 

lightest neutralinos produced by lepton flavor conserving decays of squarks and gluinos, 
as well as low-energy proton remnants. The LFV branching ratio Br(%2 — ► h^JXi) I s f° r 
example calculated in |38j in the framework of model-independent MSSM slepton mixing. 
In general, it involves a coherent summation over all intermediate slepton states. In our 
numerical calculation we use the leading order partonic cross sections a{pp — > qq), a(pp — > 
qg) and a(pp — > gg) folded with the CTEQ6M parton distribution functions together with 



(23 ) in order to calculate the number of events for the LFV process pp —> xi+X — >■ e/i+Xi+X 



expected at the LHC [391 HQ 



4 Textures 

In the following, we will consider the lepton mass matrix textures from the list of 72 types 



given in [22J, where we normalize the textures according to ( 13 ) to obtain lepton and neutrino 
masses compatible with experimental data. This list will be called our "reference list" of 
textures. Each of the 72 types of textures in the list is characterized by a set {M^, Md, Mr} 
of relevant mass matrices. We call this set of three textures a "texture set" . It is important 
to note that the textures M#, Mo, and Mr are in general not diagonal and large leptonic 
mixing can emerge from any of the matrices Me, Md or Mr. The structure of the non- 
diagonal textures is relevant for the construction of explicit high-scale theories, such as the 
Froggatt-Nielsen mechanism or discrete flavor models, generating the hierarchical pattern of 
these textures. In fact, a large number of explicit models predicting the texture sets from 
flavor symmetries have already been found for the SM [23] and also for SUSY SU(5) GUTs 

4.1 Real Textures and Relation to QLC 

As shown in [22] , each of the 72 texture sets reproduces tribimaximal neutrino mixing in the 
neutrino sector at the 3a level along with the charged lepton mass ratios^] 

e 4 : e 2 : 1, (24) 



4 We are interested in a fit compatible with minimal SU(5) but different mass spectra, e.g. realization of 
the Georgi-Jarlskog relations [41] are possible and can be implemented just as well. 



and a normal neutrino mass hierarchy of the form 



mi : m-i : m-& = e : e : 1, (25) 

where e is of the order of the Cabibbo angle e ~ 9q ~ 0.2 and 1713 = m v w 5 x 10~ 2 eV, 
which reproduces the values of the neutrino mass squared differences in 0. In the reference 
list, My) has the eigenvalues mf and Mr the eigenvalues mf (i = 1,2,3). Their ratios are 
always of the form 

mf : mf : mf = e k : e m : e n , mf : mf : mf = e p : e q : 1, (26) 

where k,m,n,p,q are non-negative integers < 2 and < p < q. The heavy Majorana 
neutrino masses are, thus, always non-degenerate. The textures in the reference list lead to 
solar and atmospheric mixing angles that are in agreement with current data at the 3<r level. 
Additionally, the reactor angle is very small and satisfies #13 < 1° which is an attractive 
parameter range for flavor models although new global fits point slightly towards a nonzero 

#13- 

The textures give rise to a large solar angle 612 ~ 33° in a way similar to QLC This 
means that the mixing angles entering the mixing matrices U x in ^ are either of the order of 
~ e n , with a positive integer n, or they are equal to 7r/4, corresponding to maximal mixing. 
The nearly tribimaximal mixing form of t/pMNS; an d in particular the observed value of the 
solar angle, is then a consequence of combining the mixing matrices U x from different lepton 



sectors, such as the charged lepton or left-handed neutrino sector, into £/pmns y i a (10)- 

Lepton mixing angles with positive integer powers of e are motivated by the observed 
CKM mixing V us ~ e, V c b — e 2 , and V u b — e 3 , and by yU-r-symmetry [12] (for more recent 
studies on /x-r-symmetry see, e.g. [13]). Differently from most applications of QLC, however, 
we do not require that Ue be of a CKM-like mixing form (i.e. that Ui ~ Vckm) or that U u 
be of the bi-maximal mixing type. Instead, in our reference list, the observed large solar 
and atmospheric mixing angles in £/pmns can originate from maximal mixing among any 
two generations of charged leptons, left-handed, or right-handed neutrinos. For instance, we 
have in 10% of the cases trimaximal mixing of the left-handed charged leptons, i.e. all three 
mixing angles in Ue are maximal. Such general forms of textures as given in the reference list 
are useful for the construction of new explicit models that can explain the observed lepton 
mass and mixing parameters. 

Consider, for example, texture set No. 1 from the reference list [22] . The charged lepton 
texture reads 

/-ae 4 

M E = ^ ae 4 e 2 -e 2 I , (27) 
while the neutrino mass matrix textures are 

M D = -^C|Ie 1 e 2 |, M R = r ^\ae 3 1 + e -1 + e I , (28) 



V2 



ae —1 ae 



be 2 


ae 3 


ae 3 


ae 3 


1 + e 


-1 + e 


ae 3 


-1 + e 


1 + e 
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where the order-one coefficients a and b, take the values a = a/2, and 6 = 2. The mass 



parameter ~ 10 2 GeV is determined by a normalization as described in (13). For this 
example, the ratios of the corresponding eigenvalues are 



mf : mf : mf = e 2 : 1 : e, mf : mf : mf = e 2 : e : 1. (29) 



The charged lepton and light neutrino mass ratios are as in (24) and (25 ). This texture set is 
an example where we have maximal mixing among the 2nd and 3rd generation of left-handed 
charged leptons, 9^ = 71/4, and right-handed neutrinos, 9^ = vr/4 (in the notation of Sec.[2j 
[22] . Note also that in it is the 2nd column that is dominant. The PMNS mixing angles 
for this texture set are given u - ' 



9 12 = - A -^ + 0{e% 23 = ^ + ^= + O(e"), 18 ~^. (30) 

Note that this texture set has the interesting property that #13 is suppressed by two powers of 
the Cabibbo angle. We have checked numerically that this is a stable feature under variation 
of the lepton Yukawa couplings and mass terms. After rotating to the basis where Mi is 
diagonal, the normalized neutrino Yukawa coupling matrix is to leading order in e given by 



e 2 _ L7e 2 _ .28e 2 \ ( R ,i (e 2 e 2 e 2N 
Y D = 0.045 I e 0.71e -0.71e ~ K m - m s) 2 £ £ £ . (31) 
-0.286 2 -1 -1 J VSm ^ V 2 1 

where we have taken m v = 5x 10~ 2 eV, m§ = 2.5 x 10 12 GeV and tan /3 = 10. With this choice 
the other heavy right-handed neutrino masses are mf = 10 11 GeV and mf = 5 x 10 11 , GeV. 



The dimensionless coefficient 0.045 multiplying the matrix in (13 X h is set by the scales m u , m R 



3 ■ 



and v sin/3. Note that even after rotating to this basis, the factors multiplying within the 
matrix the powers of e are all approximately of order one. 

4.2 Complex Textures 

The reference list contains only real matrices and therefore describes only CP-conserving 
cases. We introduce CP-violation for the different texture sets by multiplying each mass 
matrix element in the Lagrangian in ^ by an individual random phase. This means that 
we modify for the texture set {Mg, Mp, Mr} each mass matrix element to 

(M x )ij -> (M' x )ij = exp(i^) • (M x )ij (x = e, D, R), (32) 

where we assume that the 9 + 9 + 6 = 24 phases a?- vary independently on a linear scale over 
the whole interval afj G [0, 2ir[. Notice that one can rotate away 10 of the 24 phases by phase- 



redefinitions, but we will not make use of this freedom here. Moreover, the mapping in (32) 



leaves the moduli of the mass matrix elements invariant but changes the matrix R in (17) 



the PMNS mixing parameters, and the Dirac and Majorana phases 5, 0i, 02- For each of the 



5 Note that future experimental bounds might afford it to consider corrections to the solar and atmospheric 
mixing angle stemming e.g. from the charged lepton sector. 
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72 real texture sets {Me, Md, Mr} from the reference list we will thus obtain corresponding 
CP-violating texture sets {M' t , M' D , M' R }. In the following, we will call a complex texture set 
obtained in this way from the texture No. i (i — 1, . . . , 72) in the reference list a "complex 
type-i texture set" . 

Among these complex texture sets, we further consider only those with C/pmns — ^hps 
(up to Majorana phases), i.e. we study only complex cases where the PMNS matrix remains 
of nearly tribimaximal mixing form after inclusion of the CP-violating phases. Specifically, 
we require that the CP-violating textures lead to solar and atmospheric mixing angles within 
the current la errors in (J5]) while the reactor angle is very small and satisfies the constraint 
6 13 < 5° (cf. Sec. |4.l[ )J^| Moreover, we demand that the complex textures approximately 



reproduce (in our examples up to relative factors of 1.5) the mass eigenvalues in (24), (25) 



and (26), of the original texture set from the reference list they have been derived from. In 
this way, each of the 72 real examples from the reference list serves as a "parent" to a whole 
class of complex "daughter" textures that have mass ratios that are similar to those of the 
real parent texture set and also show all nearly tribimaximal lepton mixing, but they vary 
strongly in the CP-violating phases appearing in the Lagrangian. 

In addition, to ensure perturbativity of the Higgs sector, we will require that the Dirac 
neutrino matrices always satisfy 

|}f | 2 /(4tt)<0.3, (33) 

where \Y®\ 2 is the absolute value of the largest eigenvalue of Y^Yd, which is usually set 
by the largest Yukawa coupling in Yd. Note that this perturbativity constraint becomes 
particularly important for large > 10 14 GeV. 



In the parametrization (16), the complex textures will in general lead to a complex R 
matrix and, therefore, to nonzero parameters t/i. Consider, as an example, this parametriza- 
tion for a random complex texture set that is obtained from texture set No. 1 by introducing 



CP-phases via the mapping in (32). Similar to jB], we approximate M^ iag « mf diag(0, 0, 1) 
and M^ ag w m i/ diag(0,e, 1). Going to exact tribimaximal mixing £/pmns — ^hps, we thus 



obtain by comparison with ( 15 ) that 



( 771 771 ^ 1 2 

(Yuh ~ 1 v - 3 j qA, [sie^V^ H ps)i2 + ci(£/ H ps); 3 ] • (34) 
v sin p L J 

From this we conclude that Cic 2 ~ a/2 and 2siC 2 v / e/v / 3 ~ 0. Therefore, s x ~ and 
comparison with (|3l|) leads to small values xi, y% ~ 0. This implies that C\ ~ 1 and c 2 — \/2. 



In contrast to this, there are no obvious strong constraints from (34) on the possible values of 
£2,0:3, and 2/3 (apart from perturbativity arguments). We therefore roughly take an average 
value sin 22 — 2/tt. For the complex textures, we then expect from C2 — a/2 also that 
sinx 2 sinh?/2 — V%- As a result, we arrive crudely at a typical value 

y 2 ^ arcsinh (jr/V2) ~ 0(1), (35) 



6 Note that the PMNS angles 812 and 623 of the complex textures are therefore closer to the best-fit values 
(at la) than the real textures from the reference list (at 3cr). In addition, we concentrate on small values of 
013 since this limit is interesting for flavor models and is near the best-fit value. 
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Figure 1: Distribution of the parameters Xi and yi [cf. (17)] for 500 complex type-1 textures 



with a heaviest right-handed neutrino mass = 2.5 x 10 12 GeV. 



where the factor 1 / \[2 in the argument is a consequence of the approximation by tribimax- 
imal mixing. From (35), we expect in the distribution of the yi a clustering of \y2\ at values 
~ 0(1). Note that perturbativity constrains the yi not to become much larger than one. 

In Fig. [TJ we show the distribution of the parameters Xi and yi belonging to 500 complex 
type-1 textures, which all satisfy the perturbativity constraints in (33) for mf = 2.5 x 
10 12 GeV. The complex type-1 textures are all obtained from the real texture set No. 1 
in the reference list by scattering only the phases in the Lagrangian. Fig. [T] displays, as 
explained above, only the parameters Xi and yi of the complex textures which are consistent 
with nearly tribimaximal lepton mixing at the la level and reproduce, up to relative factors 
of 1.5, the lepton mass ratios given in (24), (25) and (26). 

We observe that I2/2I has, as estimated above, indeed a weak clustering at values ~ 1, 
while yi stays mostly at values \y\\ < 0.3, whereas x<± and X3 essentially vary over the whole 
interval [— 7r/2, vr/2]. In addition, 1 2/3 1 can become ~ Oil) without any clear preference for 
certain values in this range. As expected, X\ is small, with values \x\\ < tt/4. We will later 
see that the distribution of parameters allows to generate easily large Dirac and Majorana 
CP-violation phases in the low-energy theory. 

The parametrization of Yd in terms of the R matrix (17) makes it obvious that the 
inclusion of general CP- violating phases via (32) leads to an increase of the LFV rates. On 
the other hand, since the moduli of the Yukawa couplings are held fixed in this mapping one 
may expect that the increase in the LFV branching ratios still remains moderate. From (35 ), 
we can estimate that the complex textures will usually have LFV-rates that are roughly by 
a factor (ix / ' \/2) A ~ O(10) larger than for the real textures. 



5 Results for LFV Rates 

In this section, we present our results for the LFV-rates of the textures with nearly tribimax- 
imal lepton mixing. First, we consider the LFV-rates for our reference list of real textures, 
i.e. the CP-conserving case. Then, we turn to the complex textures generalizing the refer- 
ence list by including in the most general way random CP-violating phases at the Lagrangian 
level. As input parameters for the LFV-rates we take the GUT scale Mx = 2.5 x 10 16 GeV 
and an effective neutrino mass scale m v = 5 x 10 -2 GeV. Throughout, the light neutrinos 
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Figure 2: Br(/i — > cy) (top left), Br(r — > fiy) (top right), Br(r — > ej) (bottom left), and 
Br(/i — > e7)/Br(r — > fij) (bottom right) as a function of the heaviest right-handed neutrino 
mass mf for the reference list of 72 real texture sets in the SUSY scenario SPSla'. The solid 
(dashed) lines represent the current (future) experimental bounds on the respective LFV 
branching ratio. 



have the normal hierarchical spectrum given in (25). Unless stated otherwise, we will always 



refer to the mSUGRA benchmark scenario SPSla' [35]. The scenario SPSla' has a universal 
gaugino mass mi/2 = 250 GeV and a universal scalar mass m = 70 GeV at the GUT scale, 
tan/3 = 10, a positive sign sign(/i) = + of the Higgs mixing parameter fi, and a universal 
trilinear coupling parameter Ao = —300 GeV. 



5.1 Real Textures 

Let us consider first the LFV decay rates for the 72 texture sets in our reference list. All the 
textures in the reference list describe real lepton mass terms and are, thus, CP-conserving. 
Fig. [2] shows the LFV branching ratios Br(/z — > ey), Br(r — > /ry), Br(r — > cy), as well 
as the ratio Br(yU — > e7)/Br(r — > ej) for the 72 texture sets as a function of the heaviest 
right-handed Majorana neutrino mass mf. While Br(/i — > &y) varies for fixed by at 
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Figure 3: Br(/i — > erf) for all 72 real texture sets from the reference list for a fixed heaviest 
right-handed neutrino mass = 4.5 x 10 13 GeV (corresponding to Br(r — >■ /ry) « 10~ 8 for 
all texture sets) at the SUSY benchmark point SPSla'. The hierarchy in the heavy right- 
handed neutrino spectrum of the respective texture set is denoted as dark blue (strictly 
hierarchical) and light red (degeneracy between the two heaviest right-handed neutrinos). 
The symbol shape denotes the hierarchy of the Dirac neutrino mass eigenvalues: diamond 
(normal hierarchy), circle (m® is largest eigenvalue) and square (mf is largest eigenvalue). 
The solid (dashed) line represents the current (future) experimental bound. 

least one order of magnitude for different texture sets, Br(r — > /ry) changes, in comparison, 
hardly at all. As these two LFV rates scale similarly for small values of <C Mx, Br(/j — > 
Ijj) oc (mf ) 2 , in most of the textures, the ratio Br(/x — > e7)/Br(r — > /ry) is approximately 
independent of in the shown regime. Only for a handful of textures, the simple scaling is 
violated in Br(/i — > ej), and Br(/x — > e7)/Br(r — > /ry) anomalously increases with increasing 
mf. This occurs most pronouncedly in textures No. 14 and 16 which exhibit a strongly 
suppressed /t — > ef rate (cf. Fig. [3]). The r — > ej decay rate varies at least by two orders of 
magnitude, and for some textures exhibits a pronounced kink where the rate is drastically 
suppressed due to accidental cancellations. All texture sets from the reference list can have 
as large as ~ 2 x 10 13 GeV in agreement with the current bound on Br(/x — > ej) < 1CT 11 . 
Future non-observation of // — > ef at the upcoming PSI experiment, however, would imply 
for almost all texture sets that mf < 10 13 GeV. While Br(r — > fij) could be measured in the 
future for > 5 x 10 13 GeV, a measurement of r — > cy is out of reach for mf < 10 15 GeV 
for most of the textures. 

Next, let us try to gain some rough understanding of the relative size of the LFV rates in 
terms of the lepton mass hierarchies and the radiative effects of the right-handed neutrinos 
on the slepton mass matrix. In Fig. [3] we show Br(/i — > cy) for the 72 real texture sets from 
the reference list for a fixed right-handed neutrino mass = 4.5 x 10 13 GeV. It should be 
noted that this choice for is different from the value used in the rest of the paper and was 
chosen here to have Br(r — > /17) near the future experimental sensitivity. Fig. [3] demonstrates 
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Figure 4: Br(/i — > cy) and Br(r — > fi'y) for the texture set No. 1 with = 2.5 x 10 12 GeV 
in the SUSY scenario SPSla'. The green diamond shows the branching ratios for the CP- 
conserving case while the blue points denote the branching ratios for 500 random complex 
type-1 texture sets. The solid (dashed) line represents the current (future) experimental 
bound on Br(/i — > ej). 



that if Br(r — > /ij) > 10 -8 is observed in the future, all 72 texture sets can be probed at the 
SUSY benchmark point SPSla' by measuring /i — > ej. Assuming a sufficiently exact and 
independent measurement of the SUSY mass spectrum, e.g. at the LHC or by correlating 
the LFV decay rates with the anomalous dipole moment of the muon [35], the observation of 
t — > fij would provide an unambiguous determination of the heaviest right-handed neutrino 
mass m§ > 4.5 x 10 13 GeV, and an observation of fi — > ey would then constrain the number 
of viable texture sets, cf. Fig. |3j Whereas a unique identification of a particular texture set 
seems in general difficult, due to the fact that several texture sets have coincident or similar 
// — > cy rates, such an optimal observational picture would allow a much deeper insight into 
the possible leptonic flavor structure. An observable rate of r — > at future experiments 
would require a very high right-handed neutrino mass scale m§ > 4 x 10 14 GeV, which is 
incompatible with the non-observation of r — > fi'y. 



5.2 Complex Textures 

Now, let us turn to the complex case, which we will study for a specific example texture 
set. In what follows, we will focus on texture set No. 1 from the reference list discussed in 
Sec. [4], but our results are at least qualitatively representative for the complete reference list. 
Figure [4] shows Br(/i — > cy) vs. Br(r — > /ry) for 500 complex type-1 textures in the SUSY 
scenario SPSla'. The 500 complex examples are generated from the real texture (denoted 
by the green diamond in Figure [4]) in the reference list as explained in Sec. 4.2| They give 
all nearly tribimaximal lepton mixing at the la level and exhibit, up to relative factors of 
1.5, the lepton mass ratios listed in (24), (25), and (26). We observe from Fig. [4] that the 
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Figure 5: Br(/i — > cy) vs. the PMNS mixing angles and Dirac CP-phase for the complex 
textures scanned in Fig. [TJ The green diamond shows the values for the CP conserving case. 
The solid (dashed) line represents the current (future) experimental bound. 



complex textures have maximal values for Br(/i — > e'y) up to a factor ~ 200 larger than for 
the corresponding real case. In exceptional cases, Br(/i — > cy) can be smaller by three orders 
of magnitude than in the CP-conserving case. On the other hand, Br(r — > /17) only varies 
by a factor of ~ 3 due to arbitrary CP phases. Because of the potential increase or decrease 
of Br(r — > /ry) by 2-3 orders of magnitude for nonzero CP phases, it seems generally not 
possible to discriminate the complex texture sets through the observation of LFV rates. For 
completeness and consistency, we have checked that the electron electric dipole moments for 
the complex textures are several orders of magnitude below current bounds, as expected in 
general mSUGRA seesaw models with a mild right-handed neutrino mass hierarchy. |4b'j . 

Next, let us consider the LFV branching ratios as a function of the low-energy PMNS 
mixing parameters. Fig. [5] summarizes the nearly tribimaximal PMNS mixing parameters 
and Dirac phase belonging to the complex textures scanned in Figs. [T] We see that the solar 
and atmospheric angle 6*12 and 6*23 populate the complete allowed ranges 6> 12 = 33° ±1.5° and 
823 = 45° ±4° with no clear preference for any value. The same applies to the Majorana and 
Dirac CP- violation phases 0i, 02 arid S, which can take any value in the intervals 0i,2 G [0, tt[ 
and S G [0, 2n[ (shown is only the Dirac CP- violation phase 5). In particular, the complex 
textures do not prefer large or small values for the phases, even though they have been 
generated by starting out with a real texture set. Concerning the reactor angle #13, however, 
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Figure 6: Br(/i — > erf) as a function of #13 for the reference list of 72 real texture sets (left) 
and for 25 random complex type-1 texture sets (right) in the SUSY scenario SPSla'. The 
25 blue (dark) curves for texture No. 1 corresponding to a random choice of phases in the 
Yukawa couplings, while the green (light) curve corresponds to the CP-conserving case. The 
solid (dashed) line represents the current (future) experimental bound. 



we find that the presence of nonzero phases drives #13 away from very small values #13 < 1°, 
to larger values up to 5°, where we have set our selection cutoff, are preferred. This may be 
related to the way in which a small reactor angle can emerge for the textures No. 1 in the 
reference list and would therefore be the result of a possible underlying high-energy theory 
of flavor. 

It has been pointed out previously that there may be a correlation between Br(/x — > erf) 
and the value of the reactor angle 6*13 for a special choice of parameters [27J HH]. Fig. |5j 
however, does not exhibit any clear correlation between Br(/x — > erf) and any of the PMNS 
mixing angles and phases. In particular, Br(/i — > ej) varies by at least four orders of 
magnitude over the interval #13 G [0,5°]. 

Figure [6] demonstrates the lack of a universal correlation between Br(/z — > ej) and the 
value of the reactor angle #13 in more detail. The left plot shows Br(/i — > ej) as a function 
of # 13 for all 72 real texture sets, whereas the right plot shows the same dependence for 25 
random complex type-1 texture sets. The plots were generated by calculating the R matrix, 
cf. (17), and heavy right-handed neutrino masses in the respective texture (real or complex) 
and introducing a variable reactor angle through the ?7pmns matrix in (16), while keeping 
all other terms fixed. The neutrino Yukawa coupling matrix Yd is then used to calculate 
Br(/x — > erf). While the real texture sets exhibit a qualitatively similar dependence on 9 13 , 
Fig. [6] (right) clearly proves that there is no general correlation between mSUGRA LFV 
rates and the neutrino mixing angle #13 due to the unknown high-energy parameters as long 
as only the low-energy neutrino parameters are known. We have also checked other complex 
textures from the reference list with qualitatively similar results as for the complex type-1 
texture set. An analysis supporting this result but is based on a parameter scan of the R 
matrix, can be found in [49]. However, particular models can exhibit indeed a #13 dependence 
or correlation among low-energy phases [5D] . 

Lepton flavor violating rates for decays of supersymmetric particles in the context of 
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Figure 7: Contours of Br(/i — > cy) and the number of LFV events N(pp — > x° + X — > 
e/ix? + X) at the LHC (L = 100 fb _1 ) in the (m — wii/2) parameter plane. The other 
mSUGRA parameters are chosen such that Aq = —300 GeV, tan j3 = 10 and > 0. On the 
left (right) side the neutrino sector is given as in texture 1, with the CP phases chosen such 
that maximal (minimal) Br(/i — > cy) is achieved for fixed mf = 2.5 x 10 12 GeV. The dark 
(red) areas are excluded by direct SUSY searches. 



seesaw models with real parameters have already been discussed in the literature [101 EH 
l52| 153] . Therefore, we concentrate here on the effect of phases. In Fig. [7j we show for 
two specific complex examples belonging to the texture set No. 1 in the reference list the 
contours of Bi(fi — > ej) and the number of LFV events N(pp — > X2 + ^ — ^ e f 1 Xi +X) at the 
LHC [L = 100 fb _1 ) in the (m — m 1 / 2 ) parameter plane. We consider those textures that 
lead (among the complex textures in Fig. ^ to the highest and lowest rate for Br(/i — > cy). 
In the latter case, Fig. [7] (right), the MEG experiment sensitive to Br(/i — > cy) ~ 10 -13 
will only be able to find a signal in case of a rather light SUSY spectrum in the region 
r mo, m i/2 ^ 200 GeV. In the first case, Fig. Pf] (left), the part of the parameter space with 
m i/2 ^ 200 GeV is already excluded from Br(/t — > ej) < 1.2 x 10~ n . The MEG experiment 
can probe a sizable part of the parameter space which partly overlaps with the region where 
one expects a significant number of LFV decays of X2 a ^ ^ ne LHC. Note, that a subspace of 
the parameter space which is probed by the LHC is complementary to the region probed by 
the MEG experiment. 

6 Summary and Conclusions 

In this paper, we have considered the lepton flavor violating decay branching ratios Br(yU — > 
ej), Br(r — > fij), Br(r — > ej), as well as the LFV process pp — > efiXi + X in mSUGRA for 
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a broad class of lepton mass matrix textures that give nearly tribimaximal lepton mixing. 
The neutrino masses are normal hierarchical and become small due to the canonical type-I 
seesaw mechanism with non-degenerate, mildly hierarchical right-handed neutrino masses. 
The textures exhibit large leptonic mixing that can originate from the charged lepton, the 
left- or the right-handed neutrino sector. We have studied the CP-conserving and the most 
general CP-violating case obtained by varying the CP-violating phases in the Lagrangian. 
In doing so, we focused on the SUSY benchmark scenario SPSla'. We have determined 
the LFV decay rates for 72 qualitatively different sets of real texture patterns with a right- 
handed neutrino mass scale between 10 11 GeV and 10 15 GeV. These 72 textures are taken 
from a previous systematic scan of real textures generating nearly tribimaximal mixing in a 
way similar to QLC. 

We have studied in detail the LFV rates for complex textures that are obtained by 
randomly varying the CP-violating phases in the Lagrangian in the most general way. The 
resulting complex textures exhibit nearly tribimaximal PMNS mixing angles and reproduce 
realistic lepton masses with a normal hierarchical neutrino mass spectrum. All complex 
textures were, furthermore, subject to the requirement of having a perturbative Higgs sector. 
As expected for the SUSY seesaw mechanism, the electron electric dipole moment is several 
orders of magnitude below current bounds. We estimated and confirmed that promoting 
a real to a complex texture can lead to an enhancement of Br(/x — > e / ~f) of more than two 
orders of magnitude. 

For the SUSY scenario SPSla' and a heaviest right-handed neutrino mass of the order 
2.5 x 10 12 GeV, the complex textures lead to rates Br(/i — > &y) around the current experimen- 
tal bound. The rare decay r — >■ /ry could be measured in the future, whereas Br(r — > ey) 
seems to be out of reach of planned experiments. The predicted value of Br(r — > fi'-f) is 
approximately independent of the chosen texture set and the CP phases. Potentially, it can 
therefore be used to determine the right-handed neutrino mass scale quite robustly if the 
SUSY mass spectrum is sufficiently well known. We have also compared the radiative LFV 
decay // — >■ ey with the event rate for the process pp — > efix® + X via LFV decays of the 
second lightest neutralino x® a ^ the LHC. A sizable part of the mSUGRA parameter space 
that can be probed by this process at the LHC is found to be complementary to the one 
probed by the MEG experiment. 

For the complex textures, we could not find a clear correlation between Br(yU — > ey) and 
the PMNS mixing angles or phases of the low-energy neutrino sector such as the reactor angle 
#13. The solar and atmospheric mixing angle as well as the two Majorana and the Dirac CP- 
violating phases are distributed over the whole allowed range without any clear preference 
for any angle or phase. Though present, only a few complex textures give extremely small 
values of the reactor angle 9 13 < 1°. 

We thus conclude that the inclusion of random CP-violating phases at the Lagrangian 
level typically erases possible correlations between the LFV BRs and the PMNS parameters 
such as #13 that may occur for special points in parameter space. An interplay between 
a future measurement of the reactor angle and LFV BRs therefore seems highly model- 
dependent and would, even if the moduli of the Yukawa couplings are all fixed, require a 
detailed understanding of the possible CP-phases in the Lagrangian. 
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